
(12) United States Patent 

Spencer 



H 



US006356646B1 

(10) Patent No.: US 6,356,646 Bl 
(45) Date of Patent: Mar. 12, 2002 



(54) METHOD FOR CREATING THEMATIC 
MAPS USING SEGMENTATION OF 
TERNARY DIAGRAMS 

(76) Inventor: Clyde H. Spencer, 11331 McKellar 
Dr., Sonora, CA(US) 95370-9078 

( * ) Notice: Subject to any disclaimer, the term of this 
patent is extended or adjusted under 35 
U.S.C. 154(b) by 0 days. 

(21) Appl. No.: 09/252,718 

(22) Filed: Feb. 19, 1999 

(51) Intel/ GO6K9/O0 

(52) U^. CI 382/103; 348/144 

(58) Field of Search 382/103; 348/144, 

348/145, 146, 147, 149 

(56) References Cited 

U.S. PATENT DOCUMENTS 

3,829,218 A 8/1974 Alayanak 356/74 

3,984,671 A 10/1976 Fletcher et al 235/181 

4,037,048 A 7/1977 Walker 358/113 

4,139,862 A ♦ 2/1979 Haskell et al 348/3 

4,167,729 A 9/1979 Christenson et al 340A46.3 

AH 

4,678,911 A 7/1987 Sundberg et al 250/253 

4,766,551 A 8/1988 Bfigley 364/498 

4,885,697 A 12/1989 Hubner 364/497 

5,206,918 A 4/1993 Levene 382/17 

5,479,255 A 12/1995 Denny et al 356/319 

5,790,188 A * 3/1998 Sun 348/144 

OTHER PUBLICATIONS 

Anon (1996), ENVI User's Guide, The Environment for 
Visualizing Images, Ver. 2.5, Research Systems, Inc, Boul- 
der, CO, pp. 4-35 to 4-37. 

HaraUdc, Robert M., et al. (ed.), (1983), Ch. 18: Pattern 
Recognition and Classification, in Manual of Remote Sens- 
ing (Colwell, R. N., ed.), ASPRS, Falls Church, VA, pp. 
793-805. 

Chavez, Pat S. (1989), Radiometric Calibration of Landsat 
Thematic Mapper Multispectral Images, PE&RS, vol. 55, 
No. 9, Sep. 1989, pp. 1285-1294. 



King, Thomas A., et al. (1987), "Color Space Mapping 
Using Temary/Chromaticity Diagrams," in Proc. 5th The- 
matic Conf. Remote Sensing for Exploration Geol., ERIM, 
Ann Arbor, MI. 

King, Thomas A., et al. (1990), "Multispaectral Ratio Selec- 
tion Using Ternary Diagrams,*' Bull. Assoc. Eng. Geologists, 
vol. XXVII, No. 1. pp. 93-102. 

Anon (1980?), USGS MIPS documentation, http://terrawe- 
b.wr. usgs.gov. TRS/software/mips. 

* cited by examiner 

Primary Examiner — Leo Boudreau 
Assistant Examiner — Abolfazl Tabatabai 



(57) 



ABSTRACT 



A method and process is disclosed for computer processing 
and display of remote sensing, multispectral imagery for the 
purpose of identifying ground targets and classifying the 
imagery to create thematic maps. The successful application 
relies on the use of images from optimal bands; those bands 
include a near-infrared, mid-infrared, and visible band. The 
percentages of each of the three bands are plotted on a 
ternary diagram (34), which is used as a graphical device to 
allow target identification. Additionally, the ternary diagram 
feature space is segmented to define thematic classes and 
aUow thematic classification with a computer. Vegetation 
(10) and mineral clusters (11) are separated easily. Further 
discrimination within the separate vegetation (10) and min- 
eral (11) point clusters is possible. A spectral library (49), 
convolved to the bandwidth of the employed multispectral 
sensor bandpasses, and converted to ternary percentages, is 
utilized to locate specific spectral targets. The graphical 
ternary diagram (34) and false-color multispectral image 
display (46) are linked in real-time through a lookup-table 
(43) to allow an operator to interactively alter each of them 
through a change to the other. Included is a method to 
accomplish automatic, hierarchical classification of the 
images without operator intervention. The percent of veg- 
etation ground-cover (i.e. a vegetation index) can be 
estimated, and atmospheric scattering and absorption effects 
are empirically corrected, interactively. 



12 Claims, 8 Drawing Sheets 
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METHOD FOR CREATING THEMATIC characteristics of any or all of the existing clusters, exceed » 

MAPS USING SEGMENTATION OF certain thresholds, either hardcoded by the software manu- 

TERNARY DIAGRAMS facturer or selected by the operator, clusters can be split or 

lumped together. After a pre-selected number of iterations, 

BACKGROUND OF THE INVENTION 5 or alternatively a pre-selected number of clusters is attained. 

the program halts and displays the spectrally dissimilar 

1. Techmcal Field clusters, wherein the clusters are commonly pseudo-colored 
This invention relates to the computer processing and to help in their visual differentiation. The operator must then 

display of digital, multispectral imagery for the purpose of subjectively decide what every cluster represents in terms of 

identifying ground targets and classifying the imagery to lo a ground target or type of ground cover. This is again 

create thematic maps. accomplished through either photointerpretation procedures, 

2. Description of Prior Art or the collection of ground-truth from oosite inspection, or 
Ground targets in imagery have been identified commonly ^ combination of both. 

through the subjective process of an art referred to as Numerous band-ratios and indices, utilizing two bands, 

photointerpretation. This process relies on shapes, tones, have been in use for years to accentuate spectral features, 

textures, colors, and associations in the image to infer the Thresholding a band-ratio (or vegetation index) can provide 

nature of the ground feature. The results are highly corre- a binarized (two-level) classification of the images. Whereto 

lated with the skill and experience of the photointerpreter. establish the threshold is subjective. 

Multispectral imagery is a collection of coregistered It can be seen that the subjective judgement of the 

images (typically <2^), each image being collected in a operator is extremely important in obtaining accurate results 

different, broad region of the electromagnetic spectrum for all the methods discussed above. Their experience is 

commonly called a band. These bands are collected nomi- extremely important in making the subjective decisions. In 

nally between 0.400 ^on (micrometers) and 2.500 /an wave- addition, the availability of accurate ground truth is impor- 

length; this region is generally referred to as the visible/ tant in achieving good classification results. Although the 

infrared-red (VIR) region. The computer classification of algorithmic processing of the imagery can be performed for 

multi-spectral imagery, to create thematic maps, has usually an n-dimensional feamre space where n is greater than 1, the 

been of two general types, referred to as supervised and operator is constrained to only being able to see the effects 

unsupervised. of three bands at any one time. This affects the ability to 

In the case of supervised classification, regions within the either correctly select homogeneous training sites or prop- 
image, that are believed to represent relatively homogeneous erly identify the corre^onding type of ground cover for a 
clusters of pixels (picture elements) that are characteristic of particular cluster-class derived from unsupervised classifi- 
a particular type of ground cover of interest, are selected to cation. One does not know a priori what the best three 
be what are known as training sites. The selection can and multispectral bands are to help identify the various types of 
most often does proceed by an experienced photointerpreter ground cover in the image, 

marking the boundaries of the training site with an on-screen An additional problem is that one can expect that every 

cursor controlled with some pointing device such as a algorithm chosen will result in a classification of the ground 

mouse, trackball, joystick, lightpen, or keyboard arrow- cover that is different from the results of every other 

keys. The selection of partictilar training-site locations algorithm. The initial selection of an appropriate algorithm 

within an image, and their resultant shapes, is subjectively ^ from the plurality of choices available will be influenced by 

determined by an operator. The selection is based on their the experience of the operator. The acceptance or rejection 

interpretation of the image and possibly is guided by infor- of the classification results will also be a subjective decision 

mation called "ground truth*' obtained from on-site recon- made by the operator. While there are objective tests of the 

naissance on the ground. accuracy of classification, even the interpretation of the tests 

These training sites are then used by various computer 45 ^ ultimately subjective, 
algorithms, again commonly selected by an operator from a A common problem in multispectral classification, par- 
plurality of choices, and typically statistical in nature, (such ticularly with high spatial-resolution imagery, is the incor- 
as the Maximum Likelihood classifier) to sequentially test poration and ambiguous identification of a shadow class, 
the spectral characteristics of every pixel in the image Related to this is, if the operator does not recognize the 
against the pixels in the selected training sites. The tests 59 existence of one or more spectrally unique types of ground 
performed are for a similarity metric for each pixel with a cover, those classes will be relegated to an 'unknown* or 
pluraUty of multispectral bands that may represent all bands 'undefined* class with supervised classification, 
recorded or a sub-set selected by the operator. The definition Another approach for creating thematic maps, which has 
of the agreement or similarity varies with the algorithm evolved in recent years, is an extension of multispectral 
selected. The intent is to categorize or classify every pixel as 55 remote sensing. If a sufficient number of as discrete bands of 
being most like the pixels in a particular training site, when the electromagnetic spectrum are sampled, the resulting plot 
considered from the perspective of the n-dimensional feature of the apparent reflectance of a ground target approaches a 
space defined by the number of bands selected for the continuous, smooth curve. The imaging instruments that 
comparison. collect a very large number of bands, usually between about 

In the unsupervised classification approach, an operator 60 2*to 2® narrow bands, nominally between 0.400 /ma and 

selects a particular algorithm from a plurality of choices 2.500 fiin wavelength, are called hyperspectral sensors. The 

commonly provided by the manufacturer of the image advantage of hyperspectral imagery is that, with careftil 

processing software. That algorithm then performs iterative calibration, correction of illumination variations, and com- 

comparisons in n-dimensional space, where n is the nimaber pensation for atmospheric absorption and scattering effects, 

of multispectral bands chosen by the operator to be consid- 65 an apparent reflectance can be derived from the radiance 

ered The comparisons result in every pixel being assigned values for every pixel; then the pixels can be compared to a 

to a cluster of spectrally similar pixels. If the statistical library of laboratory-derived reflectance spectra for a match. 
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iThereby, ideally,..the need, for .subjective judgement on the . -Art" above, several objects and -advantages of my. invention, : 

part of the operator is eliminated. However, the atmospheric hereinafter referred to as "Ternary Multispeclral Analysis," 

corrections are not a trivial undertaking, the available spec- are as follows: 

tralUbraries are reasonably complete for minerals only, and provide a means to identify ground reflectors 

the hyperspectral sensors are not extensively deployed as 5 vvithout need for 

yet. The analysis and subsequent processing of hyperspec- ^ ^^^^^ photointerpreter; 

tral imagery to create thematic maps is generally more ./^ 

difficult and requires a higher level of skill and experience ^ * groun - , 

than for multispectral imagery (^) ^^^^ three, broad bands available on the Landsat 

Multispectral (and hyperspectral) images are invariably 10 ^W^r (TM) or similar multispectral sen- 

displayed as a single false-color image created from a , 

combination of three particular bands out of the total number (2) To provide a means to locate spectral targets defined 

of bands recorded. However, most remote sensing image as particular types of ground cover, 

processing softv>rarc currently available only provides the P) T° P'"^^"" ^ mteractively classify multi- 

abiUty to plot scattcigrams from two bands simultaneously is ^P'^'^'^''' Y^"^"^ "fj^"^ """P* °* S™""'' 

. , ,. j> J ^ . , 1 *\ • u- 1 cover without need for: 
A ternary diagram (a.k.a. tnangular plot) is a graphical 

plot based on the use of an equilateral triangle; it is used for (V a highly skilled photouiterpreter; 

displaying the relationship between three variables that are (h) ground truth; 

related in a manner such that they sum to unity or 100 (c) more than three TM bands; 

percent. For plotting purposes, the triangle is internally 20 (4) To provide a means to automatically classify multi- 
subdivided into three sets of parallel lines, each set being spectral imagery and thus gene rate a thematic map of ground 
parallel to one of the sides of the triangle. The lines, which cover without the need for: 
may be hidden for clarity, represent percentages of each of (a) an operator; 
the three components being plotted. These three variables photointerpreter; 
are commonly referred to as end-members. Each apex of the 25 ground truth; 
triangle represents 100 percent of an end-member and the / j( . t. ^ * l j 
side of the triangle opposite a particular apex represents zero (d) or more than three TM bands; 
percent. Apoint in the center of the triangle represents equal , f ^ m^;ir^ to minimize the iinpact of the 
proportions (V.) of aU three end-members. Ternary diagrams ^^;^:!' ^^^^ 'normalization; 
are tools frequently used in geology, metaUurgy, and physi- 30 (6) To provide a means to estmiate the proportion of 
cal chemistry. A familiar use is in classifying soils by end-members m a mixed-pixel; 

texture. In this appUcation, the end-members are sand, silt, P ^o provide a means to estimate turbidity and type of 

and clay. Dividing lines on the ternary diagram define ^ ^ , , . 

boundaries for soil textural names, such as "loam." Ternary .(^) P^°^^^ ^ °?^^°f ^ f^se^olor m^agery 

diagrams are exceedingly tedious to plot by hand and are 35 with constant colors for the same type of reflectors, and; 

subject to plotting error. <^^} Provide a means to interactively accomplish 

-^TTo^i'io 1-t .£CT--ijy^. . empirical corrections of atmosphenc scattcrmg and absorp- 

The U.S. Geological Survey Flagstaff Field Center wrote r j r u- i u j 

■ inort 11 Jv* -T n • tion for a multitude of multispectral bands, 

a computer program in 1980 called Mim Image Processing o^ n _4t. i_ - * j j * n u 

o . /m«VM^\ T. • -1 1.1 . .1. Li- £ ,7Av J Still further objects and advantages will become apparent 

System (MIPS). It is available to the public for VAX and , ^ r .t. • j • j 

,/ . ^ ^ , x^ino . • J 1 11 J .in from a consideration of the eosumg descnption and accom- 

Unix operatmg systems. MIPS contains a module, called ^ • • or 

"Triplot," for displaying triangular plots oif muhispcctral P > S g • 

imagery. However, there is no indication that there is any SUMMARY OF THE INVENTION 

awareness among the developers that particular band selec- ^^^^^^ ^^^^ of Ternary Multispectral Analy- 

Uons result in clusters of points that are separable into ^.^ ^^^^^^^ ^^J^ ^^^^^^^ .^^^^ 

classes of matenals. The documentation simply says, "Uses: 45 . . . f / f \ u - 

^ , . . u * J * visuahzation, image processing, mathematical techniques. 

If three multispectral images have been converted to h * i ^ * 

reflectance, then the percentage of the total reflectance ^^^^ processmg. 

represented by each can be seen in the plot." MIPS also In accordance with the present invenUon, plotting of pixel 

contains a module, called "Reflect," for converting raw data valuesfrom multispectral images from appropriate Landsat 

into ground reflectance 50 TM (Thematic Mapper) and ETM+ (Enhanced Thematic 

Research Systems Inc. (Boulder, CO) markets and dis- ^^^^^ bands (or their equivalents) on a ternary 

^ ^ \ . V , **. u r> diagram results m a feature space display that effectively 

tributes a remote sensing software package written by Better , ^ , ^ i r • i t- ... • 

Solutions Consulting LLC (Boulder. CO) called HWI. It '^•"f' vegetaUon separately from mmerals. Further, mm- 

L c * *t- * 11 «j • * 1 » *i- * 1- 1 erals are separated as a cluster along a major symmetry- axis 

has a feature, that they caU danang pixels, that links an r.u . • i j i_ i • * • . i 

image to a 2-D scatlerplot (CartesiL coordinates). Pixels » "^^"^ ^"^^ ^f*'' dispersion at approxunately 

, ^ , . • r-\ . • • 1 J neht angles to the general elongation of the mineral cluster, 

selected within a region of interest in the image arc plotted „ • i • j r . j i 

. . t . J 1 • * AdditionaUy, various kinds of vegetation are dispersed along 

on an orthogonal-axis scaUerplot and, conversely, points j • i- i * * j xi. • i 

1 ^ J ■ 1 . .t- • J - • 1 - a second maior Imear cluster, separated from the mmeral 

selected m the scatterplot cause their corresponding pixels in cluster 

the image to be highlighted. Again, however, there seems to . * . ^ ^ . . 

be no awareness on the part of the developers that one can 6" According to one feature of the invention, a ternary 

have any a priori knowledge about the identity association of ^"^^^^ '^^"^'^ from spearal reflectance libraries coo- 

thepointsiQthescatterplot.orawarenessoftbevalueofany volved to the baodwidths of particular bands froin an 

particular band combinaUons. ^^^f}"^, multispectral seizor system (e.g Undsat TM), 

will display point clusters definmg separable regions of a 

OBJECTS AND ADVANTAGES 55 plurality of landcover themes. 

Accordingly, besides the objects and advantages of mul- According to another feature of the invention, a computer- 

tispectral classification as described in "Description of Prior graphics ternary diagram has a real-time link to a computer- 
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. displayed image, through a lookup . lablcu This allows , o.^ . FIG. 7C shows vegetation and . mineral -alignment bands ^ 
selected pixels in the image to be identified as particular for TM bands 4, 5, and 3. 
points in the ternary diagram and, conversely, pixels corre- 
sponding to selected points in the ternary diagram can be REFERENCE ^am4ERALS IN DRAWINGS 
highlighted in the image. This reciprocal functionality 5 vegetation cluster 
allows the identification of targets in the image and/or the n mineral cluster 
generation of spectrally-derived thematic maps of ground- 12 bisector 
^over. 13 temperate-region plants 

According to another feature of the invention, schemas 14 arid-region plants 

can be constructed to automatically classify an image into 10 \^ desert soils 

themes. The schemas can be customized for particular goals. 16 senescent plants and litter 

In still another feature of the invention, ternary diagrams 17 non-silicate minerals 

can be used to empirically correct for variations in solar 18 felsic, rock-forming minerals 

illumination and atmospheric absorption and scattering. enstatite pyroxenes 

Other objects, features, and advantages will be readily ^0 evaporite minerals 

apparent. 21 snow and ice 

22 pure water 

BRIEF DESCRIPTION OF THE DRAWINGS 23 image data files 

wLj - t.^ 24 computer (CPU) 

In the drawmgs, cbsely related figures have the same 20 25 computer input devices 

number but different alphabetic suffixes. 26 band choices 

The various advantages of Ternary Multispectral Analysis 27 image/diagram control choice 

will become apparent to those skilled in the art after studying 28 read image data 

the following specification and by reference to the drawings, 29 test for data EOF 

in which: 25 30 choice for raw/corrected pixel 

FIG. LA shows a triangular plot of minerals and vegeta- 31 reflectance correction 

tion obtained from laboratory spectral libraries convolved to 32 calculate percentages 

the bandwidths of Landsat TM bands 4, 7, and 1. 33 plot points 

FIG. IB shows a triangular plot of minerals and vegeta- 34 display in diagram window 

tion obtained from laboratory spectral reflectance hbraries 35 monitor (display device) 

convolved to the bandwidths of Landsat TM bands 4, 7, and 36 release control 

2. 37 select diagram/library option 

FIG. IC shows a triangular plot of minerals and vegeta- ^8 select region of interest 

tion obtained from laboratory spectral reflectance libraries ^^^^ image data 

convolved to the bandwidths of Landsat TM bands 4, 7, and 40 test for data EOF 

3 41 calculate reflectance 

FIG. 2A shows a triangular plot of minerals and vegeta- calculate ^rcentages 

tion obtained from laboratory spectral reflectance libraries crea e 

convolved to the bandwidths of Landsat TM bands 4, 5, and ^ 44 test tor tolerance 

^ ' ' 40 45 highlight pixels 

1 . t t . 1 1 46 display in image window 

FIG 2B sho^vs a tnangular plot of miMrals and vegeta- ^j^^, ^^^^^^ 

tion obtained from laboratory spectral reflectance hbraries ^^^^ Ubrarv 

convolved to the bandwidths of Landsat TM bands 4, 5, and spectral library 

^- 45 50 test if all spectra read 

FIG. 2C shows a triangular plot of minerals and vegeta- 51 mineral alignment band 

tion obtained from laboratory spectral reflectance hbraries 52 vegetation aUgnment band 

convolved to the bandwidths of Landsat TM bands 4, 5, and 53 combined alignment band 

^- 54 vegetation alignment band 

FIG. 3 shows a ternary diagram with general segment 50 55 mineral alignment band 
boundaries between major reflectance targets such as 

vegetation, minerals, and soils for Landsat TM bands 4, 7, DESCRIPTION OF THE PREFERRED 

and 1. EMBODIMENT 

FIG. 4A shows a triangular plot of a typical Landsat scene ^' Overview 

for TM bands 4, 7, and 2, using raw digital numbers. ss ConventionaUy, a ternary diagram is displayed with one 

FIG. 4B shows a triangular plot of a typical Landsat scene 0/ ^^^^^ P^^^^f »° ^^^^^'^ the display medium, as 

for TM bands 4, 7, and 2, using values converted to ^"^^^^ ^ Uow^vcu the 

exoatmospheric reflectance. orientation of the triangle or the order of assignment of 

, ^. end-members only changes the appearance m the same 

FIG. 5 shows a schematic flow diagram of the processes ^^anner as rotating a photograph or observing a photograph 

mvolved m the general computer processmg. ^ ^ fundamental relationships between the data 

FIG. 6 shows a flow diagram of an automatic classifica- clusters are only rotated or reflected, not transformed, 

tion schema. Therefore, the preferred embodiment will be presented in the 

FIG. 7A shows vegetation and mineral ahgnment bands context of the conventional display, 

for TM bands 4, 7, and 5. 55 Three multispectral images, hereafter referred to as bands, 

FIG. 7B shows combined alignment bands for TM bands are selected from a plurahty of bands available from an 

2, 7, and L imaging muUispectral sensor such as the Landsat Thematic 
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Mapper (TM). Each of the pixels (pictuce.,elements),.from 
each of the three digital image bands are read (ingested) 
from the storage media by a computer under programmatic 
control. A new floating-point value is calculated for each 
pixel representing its relative percentage compared to the 5 
whole (i.e. EQ. 1: Percentage^(x,y)=DNj(x,y)/(DNi(x,y)+ 
DN2(x,y)+DN3(x,y)); Percentage2(x,y)=DN2(x,y)/(DN,(x, 
y)+DN2(x,y)+DN3(x,y)); PerccDtage3(x,y)=DN3(x,y)/(DNi 
(x,y)+DN2(x,y)+DN3(x,y)). DN„ represents the digital 
number of a pixel, usually, but not necessarily, a byte value; 
n represents one of the three selected bands, and; (x,y) 
indicates a pixel at a particular spatial location {x,y}. If the 
image data have been preconverted to reflectance, then the 
data will be floating-point values between 0 and 1, instead 
of binary bytes. 15 

It is a consequence of the characteristics of ternary 
diagrams that if two of the percentage values are calculated, 
the third is automatically determined by plotting the two 
known values. Hence, it is not actually necessary to calcu- 
late all three percentage values as above (i.e. EQ. 2: 20 
Percentage3(x,y)-100-(Percentage3(x,y)+Percentage2(x,y)) 
) to plot a point. 

A look-up table (LUT) is created where the new, deriva- 
tive ternary-percentage values (PerceotageJ for the plurality 
of pixels converted are linked to the positions of the original ^5 
pixels. The sequential, indexed position in the LUT corre- 
sponds to the relative position of the image pixel in the 
sequential image file. One may choose to display all pixels, 
all pixels from a delineated region of interest (ROI), or only 
a representative sample from either. 3Q 

Tliis is important: Not only is a ternary diagram created 
showing the distribution of pixels in the triangular feature 
space, but every point in the ternary diagram corresponds to 
a known location in the original image. Hence, one may 
either highlight pixels in the image that correspond to a 35 
particular point or region(s) in the ternary diagram or, 
conversely, select a pixel or ROI in the image and highlight 
the corresponding points in the ternary diagram. This latter 
process allows one to identify what the materials are in a 
region of interest in the image. ^ 
2. Physical Basis of the Invention 

The key to successful classification of multispectral 
remote sensing imagery is to select a feature space that 
results in distinct, separable clusters of pixels that are 
representative of different classes or themes. Particular com- 45 
binations of Landsat TM bands, plotted pixel by pixel on a 
ternary diagram, produce a strong discrimination between 
organic and non-organic reflectors. Furthermore, with these 
particular combinations, the plotted points representing min- 
eral reflectors cluster along one line bisecting the triangle, 
allowing various mineral anion and cation groups (e.g. 
silicates or borates) to be identified along the cluster. The 
plant species similarly plot in a diflftise cluster along one 
edge of the ternary diagram, providing discrimination of 
different plants such as grasses, shrubs, and deciduous or 55 
evergreen trees. The band choices are critical. Using the six 
reflective (non-thermal) bands of Landsat TM, there are 20 
different combinations (order being significant) of three 
bands that can potentially be used. However, only six 
combinations, as shown in FIGS. 1 and 2, provide useable 
plots. Other band combinations result in undifferentiable 
clusters of points, with vegetation plotting amidst the min- 
eral point-cloud. 

The optimum selection of bands to use is as follows: 

1) a near-infrared (NIR) band (e.g.TM 4); 55 

2) a mid-infrared (MIR) band (e.g.TM 5 or 7), and; 

3) a visible band (e.g.TM 1. 2, or 3). 
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> :The NIR band is principally^ used to discriminate vegeta-. 
tion types and separate them from minerals. The MIR band 
is used to discriminate minerals. The visible band is used to 
discriminate water, and to disperse the vegetation and min- 
eral clusters at approximately right angles to their general 
trend. 

The following table lists the nominal bandpasses for the 
reflective TM bands: 



Band 1 


Q.45-0.52 fan 


(blue-green) 


Band 2 


0.52-0.60 fan 


(green) 


Band 3 


0.63-0.69 pan 


(red) 


Band 4 


0.76-0.90 fan 


(NIR) 


Bands 


1.55-1.75 fan 


(MIR) 


Band 7 


Z08-2.35 fan 


(MIR) 



Bands from other sensor systems that have approximately 
the same spectral bandpass-centers as the Thematic Mapper 
bands may be used similarly. Hyperspectral sensors may 
have their bands convolved to correspond to the TM band- 
passes or a single representative, narrow band within each of 
the TM bandpasses may be used in lieu of the TM bands. If 
using hyperspectral imagery, the best discrimination for a 
particular mineral would result if a band were selected from 
within the MIR regions that corresponds to a distinctive 
absorption or reflectance feature in that mineral *s spectrum, 
if one exists. 

The French SPOT Image satellite, SPOT 4, does not have 
bands equivalent to TM band 1 or TM band 7. Therefore, the 
use of SPOT 4 XS imagery narrows the available choices to 
TM 4, 5, and 3 or 2 equivalents. The older SPOT satellites 
also do not have the equivalent of TM band 5; so, imagery 
from SPOT 1, 2, and 3 can not be used. 

For selecting the optimum triplet of bands for Ternary 
Multispectral Analysis, a nimiber of factors come into play, 
generaUy resulting in compromises. It is desirable to maxi- 
mize the band variances, while minimizing the interband 
correlations. The optimum bands to accomplish that are 
slightly different for vegetation than for minerals. Therefore, 
the nature of the classification task or the area selected will 
influence the preferred choice of bands, 
3. Detailed DesaiptioD of FIGS. 1 to 3 

FIG. lA is a ternary plot derived from various laboratory 
reflectance spectra that have been convolved to the band- 
widths equivalent to TM bands 4, 7, and 1. The vegetation 
points 10 cluster along the left edge of the triangle. The 
mineral points U generaUy cluster from the center of the 
base of the triangle upward to the apex opposite. Water can 
be expected to plot to the right of the mineral cluster, with 
pure water near the right-hand comer of the triangle. TM 
bands 4, 7, and 1 give the greatest dispersion for the mineral 
clusters, hence potentially the greatest discrimination of the 
minerals. However, there is also the greatest overlap 
between vegetation and minerals, leading to confusion 
between them. This choice of bands is the theoretical 
optimum under ideal circumstances. However, subsequent 
description wiU show why other band combinations may be 
necessary or preferred under various circumstances. 

FIG. IB is a ternary plot derived from the same laboratory 
reflectance spectra as above, convolved to the bandwidth 
equivalents of TM bands 4, 7, and 2; 

FIG. IC is a ternary plot derived from the same laboratory 
reflectance spectra as above convolved to the bandwidth 
equivalents of TM bands 4, 7, and 3. The clusters (10 and 11) 
are similar to those obtained in FIG. lA. TM bands 4, 7, and 
3 give reasonably good separation between plants and 
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minerals;iJowever,'mineral points 11 .are essentially all in a . 
fat line with only slightly greater lateral dispersion than what 
a TM 4, 5, and 3 combination provides. 

The choice of TM bands 4, 7, and 2 seems a good 
compromise between 4, 7, 1 and 4, 7, 3. Another advantage 5 
of the TM 4 and 7 combinations is that images composited 
from the ratios will show vegetation greener than 4 and 5 
bands, which will be more yellow. 

FIG. 2Ais a ternary plot derived from various laboratory 
reflectance spectra convolved to the bandwidth equivalents lO 
of TM bands 4, 5, and 1; FIG. 2B is a ternary plot derived 
from the same laboratory reflectance spectra as above con- 
volved to the bandwidth equivalents of TM bands 4, 5, and 
2. 

FIG. 2C is a ternary plot derived from the same laboratory 15 
reflectance spectra as above convolved to the equivalent of 
TM bands 4, 5, and 3. This particular band combination is 
unique because the vegetation 10 and mineral 11 clusters are 
the most Unearized. The use of TM bands 4, 5, and 3 for 
image interpretation offers the advantage that any points that 20 
lie between the vegetation cluster 10 and the mineral cluster 
11 are most probably mixed-pixels. 

The use of TM bands 4, 5, and 2 has greater mineral 
dispersion than 4, 5, 3. However, there is more overlap with 
the mineral cluster and vegetation cluster with 4, 5, 2 than 25 
with 4, 5, 3. TM bands 4, 5, and 3 also provide less 
dispersion of the vegetation than 4, 5, 2. However, TM bands 
4, 5, and 3 do have the additional unique property of 
distinctly linearizing the enstatite (hypersthene-bronzite) 
pyroxene mineral solid-solution series the best of all. 30 

The use of TM band 7 produces a longer theoretical 
mineral-cluster U than TM band 5 and moves the vegetation 
cluster 10 farther away. However, examination of 2D spec- 
tral graphs indicates that many times the relative positions of 
spectral lines are interchanged between bands 5 and 7. The 35 
mineral enstatite-offshoot line is best expressed using TM 
band 5 rather than band 7. Serpentine minerals (antigorite, 
lizardite, and chrysotile) show considerable dispersion along 
the TM band 5 axis. This offers the potential for 
mineralogical/chemical differentiation in what is usually 40 
rendered in images as a dark, nondescript lithologic unit. If 
masking is employed, so that only pixels from the serpen- 
tinite body are plotted, then it may be possible to classify the 
body into spectrally distinct regions. 

Because the two main clusters — ^vegetation 10 and min- 45 
cral 11 — are each elongated along the direction of the 
MIR-band axis, each of these clusters can be further subdi- 
vided into groups. 

FIG. 3 shows a ternary diagram plot derived from various 
laboratory reflectance spectra that have been convolved to 50 
the bandwidths equivalent to TM bands 4, 7, and 1. In this 
plot, further discrimination of the various plant and mineral 
groups is shown, along with approximate segment bound- 
aries (shown in magenta) to create a sample of possible 
thematic classes. This is not intended as the only possibility, 55 
only one example of a plurality of possibilities, determined 
by the classification requirements. 

Within the vegetation cluster 10, pine trees and other 
evergreens lie at the end of the cluster with the smallest MIR 
band values and largest ^^R values. Broad-leaf deciduous 60 
trees occupy the center of the cluster, with leafy plants such 
as commercial vegetables generally being farther along the 
line of the cluster; these are shown in FIG. 3 as temperate - 
region plants 13. Arid-region plants 14 tend to plot closer to 
the mineral cluster. Senescent plants, bark, and litter 16 can 65 
be found merging into the mineral cluster at the extreme low 
NIR, high MIR end. 
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, .Desert soils 15 are shown intermixed. with the- top.half oft), 
the mineral cluster 11 and with arid-region plants 14 and 
litter 16. Desert soils 15 plot in an elliptical cluster, parallel 
to the main mineral cluster 11, but to the left of the bisector 
12 (yellow line) and sUghtly higher than the felsic mineral 
cluster 18. 

Virtually all common rock-forming silicates (as well as 
hmestone and dolomite) cluster near the center of the ternary 
diagram. The abundant rock-forming minerals 8 such as 
feldspars, quartz, calcite, and dolomite tend to cluster very 
tightly between 30 and 40 percent for a TM band 4, 5, 3 
triad. ATM band 4,7,2 triad is only slighdye more dispersed. 
The composite plot of all felsic minerals 8 is largely 
restricted to the range of 25 to 45 percent MIR with the use 
of TM band 5. 

Most silicates plot from 20 to 70 percent along the mineral 
line for the use of TM 7 and 20 to 60 percent with the use 
of TM 5. The tectosilicates are restricted to a range of about 
20 to 40 percent MIR contribution, along the mineral line. 
However, the other siHcate structiu-al groups overly them. 

Iron-bearing silicates, sulfates 7, and oxides 7 (listed 
under "Reference Numerals'* as "non-silicate minerals" tend 
to be at the end of cluster 11 with high MIR values. 

There is a notable deviation from the general mineral line; 
it is a separate, distinct line at an acute angle to the bisector. 
This line represents the enstatite (hypersthene-bronzite) 
pyroxene solid-solution series 19. Fortuitously, it is on the 
side of mineral cluster 11 opposite from vegetation cluster 
10. It does, however, intrude into the region characterized by 
the identification of water bodies. The rest of the inosihcates 
plot upward along bisector 12. 

As a general rule, for the preferred band combinations, 
(e.g. TM 4, 7, & 2 or TM 4, 5, & 3) minerals that are 
green(i5h) to the himian eye tend to plot on and to the right 
of the bisector 12, while red(ish) minerals tend to plot to the 
left. Minerals that are white or colorless are exclusively 
below the 50 percent line of the TM 7 or TM 5 percentage; 
most are below the 40 percent Hne. 

Within mineral cluster 11, evaporite minerals 20 such as 
borates and hahdes occupy the end of the cluster with the 
smallest MIR-band values. Borates plot from zero to 20 
percent along the mineral-cluster line for the use of TM 7 
and 5 percent to 30 percent along the line for the use of TM 
5. These values represent the percentage contribution of the 
MIR band. Some carbonates also plot between 5 and 15 
percent along the line with TM 7, ahhough most cluster 
around 30 percent. Halidcs plot from 10 to 35 percent along 
the mineral line with TM 7 and from 20 to 35 percent with 
TM 5. Sulfates are restricted to plotting between 0 and 42 
percent along the mineral line for TM 7 use; they extend 
from 5 to 50 percent along the mineral line for TM 5 use. 
Generally, points that plot at the base of the mineral line (low 
MIR values) can be expected to be evaporites 20. 

All oxides (except water and ice) plot high (20 to 90%) 
and along, or significantly to the left of bisector 12. Pure 
water 22 plots at the extreme right-hand tip of the triangle, 
represented by the percentage of the visible band. Water with 
sediment loads will plot intermediate between pure water 22 
and the mineral cluster 11. The position will be determined 
by the type and amount of sediment. Snow and ice 21 plots 
adjacent to and to the right of the evaporite minerals 20. 

Typical secondary gossan-minerals plot from about 38 to 
59 percent MIR with the use of TM 5. Sulfides have a greater 
variance along the mineral line, ranging from just imder 20 
percent to just under 60 percent with TM band 7. The use of 
TM band 5 restricts the lower range to about 29 percent. 
Noteworthy, is that the distribution of sulfides about bisector 
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, 12 (mineral line) of the triangle is approximately>symmetri- ... ..are overlapping ve gelation,- the TM.4,>7,. lb and combination 

cal for a TM band 5 plot. While for a TM band 7 plot, there is best for identifying or locating minerals. Actually, the 

is a greater dispersion on the left side towards an increasing same can be said of vegetation. Although, all of the com- 

TM band 4 percentage. binations using TM band 7 have similar diversions for 

Copper-bearing minerals are an anomalous group in that 5 vegetation, 

many of them plot high and to the right of bisector 12. TM bands 1 and 2 have more severe Rayleigh scattering 

One can expect that since monomineralic outcrops are than TM band 3 and therefore may be more difiBcult to 

rare, most mineral pixels will plot near the center of the correct in scenes. The lateral variance in the mineral points 

mineral cluster. That is because even if a uniform rock increases as one progresses from TM band 3 to TM band 1. 

outcrop were to dominate a pixel, the pixel is still inherently lO Depending on where the scene is located and what mineral 

a mixture of the minerals constituting the rock. This is a groups one is interested in, the scattering may or may not be 

problem confronting even hyperspectral systems, where a problem. However, the use of TM band 1 or 2 will 

destructive interference of the reflection spectra tend to probably do a better job of discriminating and identifying 

flatten the resultant spectrum, obscuring absorption features. water than would TM band 3. 

However, it is the uncommon monomineralic outcrops that 15 The shorter the wavelength used for the visible band, the 

are often of economic value. Additionally, the extreme width more Rayleigh scattering present, reducing the contrast and 

of the theoretical mineral cluster 11 is largely the result of thus variance. The potential problems with effective atmo- 

the influence of relatively imcommon minerals that rarely, if spheric correaion of TM band 1, and inevitable reduction in 

ever, occur as substantial monomineralic outcrops. variance and signal-to-noise ratio (SNR) may mitigate 

Since many plants and plant communities have open 20 against the practical use of TM 1. 

canopies, the soil below is often visible, resulting in pixels There is a similar penalty to be paid at the longest 

that have values that are a composite of the vegetation and wavelengths. The description has centered on the theoretical 

soil. These are commonly refenred to as mixed-pixels. The clusters derived from laboratory-derived spectral libraries, 

result is that in real-world situations, the vegetation tends to There is less radiant energy emitted by the Sun at TM band 

plot closer to the mineral cluster than one would find with 25 7 wavelengths than at TM 5, resulting in less recorded 

using laboratory spectra of the plant leaves only. variance and a consequent lower SNR that may make the 

The overlap that occurs between vegetation cluster 10 and practical application less effective than theoretical consid- 

mineral cluster 11 points derived firom laboratory spectral erations would suggest. 

libraries is typically attributable to senescent plants and 4. Example of Actual Use of Invention — FIG. 4 
litter, and dead or dying non-arid region plants 16. The 30 FIG. 4A is a triangular plot of actual raw TM data from 
region between vegetation cluster 10 and mineral cluster 11 bands 4, 7, and 2. The data are from a typical U.S. south- 
is occupied also by arid-region plants 14, soils 15, and western desert scene containing some irrigated agriculture, 
potentially mixed-pixels composed of both plants and bare The percentage marks and grid lines have been omitted. It 
rock or soil exposures. If the area being studied is not arid, displays the characteristic vegetation cluster on the lower- 
then obviously one is not looking at arid-region plants. 35 left and mineral cluster (predominantly mineral soils) along 
Therefore, the points most likely represent mixed -pixels. the vertical, with a string of connecting mixed pixels. 
One may use the position to estimate the percentage of plant Additionally, there are points in the right center that indicate 
cover. muddy irrigation water. No atmospheric corrections have 

Arid-region vegetation rarely covers more than 20% to been applied. 

30% of the ground. Therefore, the vegetation's effect on the 40 FIG. 4B is a triangular plot of TM data from bands 4, 7, 

position of soil or mineral points will be substantially less and 2, that have been converted to exoatmospheric reflec- 

than vegetation in humid areas. The mixed -pixel effect is tance for the same scene as FIG. 4A. The correction used the 

less. The mineral points are translated a shorter distance and metadata and algorithms provided by EOSAT for exoatmo- 

therefore will be easier to correctly identify. Further, the fact spheric correction- A correction of -8 DNs was applied to 

that typical desert vegetation is located intermediate 45 the TM band 2 data to correct for Rayleigh scattering in the 

between the mineral cluster and humid -climate vegetation green region of the electromagnetic spectnmi. Hie point 

give it less weight in translating soil and mineral points. clusters are very similar to FIG. 4A. 

For locating any particular mineral in an image (other 5. Description of Computer Flow-Diagram — ^FIG. 5 

than the hypersthene/bronzite species), a band combination FIG. 5 shows a general schematic flow-diagram of the 

such as TM 4, 7, and 2, which disperses the mineral cluster 50 essential elements of the disclosed methods and processes, 

into a broader cloud, is preferred over a TM 4, 5, 3 as implemented for computer use. Multispectral Image Data 

combination because there should be fewer false-positives. 23, commonly stored as files on removable media, are 

TM bands 4, 7, and 1 would be even better, but Rayleigh imported into a Computer (CPU) 24. Image Data 23 are 

scattering is most severe in the blue region of the electro- ingested by an image-processing program employing Ter- 

magnetic spectrum (TM band 1) and therefore may require 55 nary Multispectral Analysis, and stored as an array of data, 

corrections that reduce the reliability of the segmentation of An operator, using typical Input Devices 25, such as a 

the image's ternary diagram. TM band 2 is a good compro- keyboard and mouse, invokes the Ternary Multispectral 

misc for the minimal-correction band 3, (which generates a Analysis module and Selects 26 bands to be used for ternary 

thinner cloud, approaching a line in shape) and band 1 diagram display. The operator then decides whether to Select 

(which gives maximum dispersion but is most in need of 60 27 image or diagram Interaction priority, 

scattering correction). However, using either band 1 or 2 If the image option is chosen, the stored multispectral data 

reduces the certainty of apparent mixed-pixel points lying array is Read 28 one pixel at a time. At least one set of a 

intermediate between the vegetation cluster and the bisector three -band triad, as previously chosen 26, is read from the 

of the ternary diagram. They may actually be mineral array. An End Of File (EOF) Test 29 or similar check to 

outcrops. 65 verify completion of data reading is performed. Next, the 

Therefore, assuming that if one can perform good atmo- operator may make a Decision 30 to use raw data or convert 

spheric corrections, and one can accept false-positives that the band values for each pixel to apparent reflectance. This 
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-.vj:DecisioD - 30 step, is nol.-.a. ;critical . elemenU since ^Ternary .v: . - be; desirable . in - the .:same^.implementalion; Tbcr spectra are 
Multispectral Analysis can be used in this mode with either then Read 48 frona spectral Library 49. The spectra could be 
raw data or apparent reflectance data. However, the flex- entered as reflectance values, but it would save computa- 
ibilily of being able to make the choice and not be limited tional time if they were preconverted to percentages for all 
by the particular implementation adds value through 5 expected triads of bands to be used. Therefore, the 
increased functionality If it is decided to work with apparent conversion step is omitted from the flow chart and it is 
reflectance values, the next step is to Convert 31 the raw assumed that preconversion has been performed. A Check 
DNs for each pixel to apparent reflectance. Whether the last ^0 is performed to be sure that all requested spectra are 
step is omitted or not, the next process is to Convert 32 the r^tneved. At this point, starting with the Read 39 unage 
. , , , . * ■ *u 1 • -1 « array function, the logic flow follows the same path as for 
pixel values to ternary percentages in the manner explamed lO , -J^^, ^ , ' ^ . « ^ ^ 
u /rr\ i\ • .u J • *• J* r» • * the KOI OptioH 38, except that apparent reflectance Con- 
above (EQ. 1) in the description of ternary diagrams. Points . \ ^ MA * 

are then Plotted 33 on the ternary diagram Display 34 in the ^"f"" ^ J^'^^'i^^H 44 is 

graphical user interface Window displayed on a computer Perf°™ed against the Lbrary spectra selected, rather than 

Monitor 35. An alternative embodiment might be to use a segment region or point(s) fronn the ternary diagram, 

second monitor or to sequentially alternate between the 15 ^"^'/^^ of Automatic ClassiflcaUoi^FIG 6 

image and graphical display. When all of the image pixels ^ \f'°P^^ ["^^^ flow^agrain of one of a 

have been read and plotted the function will terminate and """l'''"''.* algorithmic schema that can be developed for 

Release 36 control to aUow other actions by the operator. "'PlP'ff dassjfication of multispectral miagery 

It is assumed that the above procedures have been per- with Ternary Multispecttal Analysis The schema is based on 

formed at least once so that there is a ternary diagram 20 f^"!"^/^ ,- " comparative (GRE/«ER THAN. 

Display 34 ofeitberaU the image pixels, or a sub^t selected ^ THAN) logic tests used with Boolean AND tests to 

by the operator. TTie plurality of pixels chosen is not critical ^S-^^f tn^n^\^ feature space a pnori. using ternary 

as any image may be considered a superset of smaller tiles ^'f'^'^^f- f 

or irregular polygons. However, it do« add flexibility and ^."^ '^f ™4^°=:»5 AND TM1%<21, be 

utility if an operator can select a region of interest (ROI) 25 "^te^preted as foUows: If a p«el value for the ™ band 4 

interactively f?om a displayed image and only plot points on "^^f . "'.rAA , ,^ .\ .tT^ 

the ternary diagram that a.rrespond to those pixels that Ue P^^'l" ^ If^Sf f P«.^^°'' f 

within the confines of the ROI. It allows identification of P"'^"'^ .^f ; T^! ^^'"l T ''P'"^"'"" 

. , 1 • 1 . JT 11 -.L' i_ J r *L tive of a plant and is tested further in the same manner 

those pixels exclusively that fall withm the bounds of the ? ^ ^ , v "T • i. 

j^Qj f / (downward m the nght-hand column) to determme the 

The operator then next makes a Choice 37 of a library or 8!=°"''' '''l^f ^"'^^ ' ^^^^ -T ^'f^'^'^ ' T^'' 
diagram option. If the diagram option ischosen, the operator P"f ' "T. ^^J^ representative of a mm- 
then would use cursor ktys, mouse, or some otheHnput ""^ '^en ^e tested ftirther (downward m he left-hand 
Device 25 to select one or more pixels, referred to here as ^°''^'»°) the mmeral class. AU subsequent 
a Segment 38 from displayed ternary diagram 34. If a very 35 diamond-shaped boxes work in the same manner, just using 
II L r • * • 1. * u -J • ui / different values and bands. The approach reqmres the coo- 
small number of points is dioscn, it may be desirable to ■ ^-.^ ■ . . ^ . i ^ r 

,1 Vifi* * ^ 4UTn- version of the imagery to apparent reflectance values before 

assign a tolerance level for later acceptance of a match. The .... i- j j . j i 

imaie data array will again be Read 39. Once more, EOF *f apphed. TJe order, number of tests and value 

Test 40 wiU be applied to conclude data reading and Release '^'^f 0"sl»PS are variable, bemg determmed by the classifi- 

36 control. The process of Calculating 41 an apparent 40 "''"^ ^equuements of a particular task. Tte unportant 

reflectance might be the next step, but can and should be ^'T,"".?-! '^^''T . ' hierarchical classification 

omitted if ternary display 34 is derived from raw (DN) pixel facilitates rapid compuution because it only requu^s 

values. The next step is to Calculate 42 ternary percentages ««npansons of one or two values at each step to define the 

foreverytriadseIected.Alook-uptable(LUT)isthenFilled feature-space segment. However, there is no fundamental 

^ . J i_ ^ 1 reason why the third value could not be used, e.e. to trade 

43 so that there is a correspondence between the sequential 45 *" ^ """"^ " """^ 

position in the LUT and the sequential position of the pixels ott^eed tor accuracy 

ill the image data array. The data in the LUT consists of the ^ ^"T. T"" Co^ections-FIG. ^ ^ ^ ^ 

ternary percentages for one or mote triads of bands. Tlie next , ^1°- °^ » f "^^^ '^'^^""^ bands4. 7, and 

stepistoCheck44toseeifthepixeltemarypercentagesare 5 with both mmeral and vegetation libraries plotted. A 

•*u- « J- nriT * * I c ^« mmeral alignment band 51 and a vegetation alignment band 

within the ternary diagram ROI, or acceptance tolerance for 50 52 are to used for verif in ro r atmos heric correc 

a small number of pixels. If the pixel ternary values pass the tion^ ^ proper a mosp enc correc- 

acceptance test, then the corresponding pixel in the image is r^,)-, «^ , t- prr^,. 1 ^ - 1 

Highlighted 45 in some way. One way would be to replace , "P J^u^.^"^-* diagram of TM bands 2, 7, and 

the co^ponding pixel in image Display 46 with a brighdy ^ both hbranes agam plotted Smoe the vegetation and 

colored pixel, or blink the pixel on and off, or both. This 55 "I"""' ^'"f' "f "'^b other, there is only a single 

r J* *i- TirT- alignment band 53 tor this combmation. 

iterative process of readmg, converting, filling the LUT, „ . v: , . , ^ ^ 

. ^- , 1 • • r *• *T 11 .u ■ rlG, 7C shows a ternary diagram of TM bands 4, 5, and 

testmg, and replacmg pixels continues untu all the image . , ^ . , ^ . . .1 \ 

, 1 . . 1 3 with both libranes plotted. A vegetation augnment band 54 

pixels have been examined, . . . , ^ %u m • mr- a -^^i i- 

wc ^ ^- Tu /^u • *L l u L J (coincident with 10 m FIG. 2C) and a mmeral augnment 

If at diagram -library Choice 37, the library option had > . .t^^.;-.T^-.^^ . 

. u .u * ij u J * c* I * band 55 (comcident with U m FIG. 2C) are shown, 

been chosen, the operator would be expected to Select 47 60 ^ ^ 

one or more entries, by name, from a Library 49 of spectra. PREFERRED EMBODIMENT— OPERATION 

These selections could be mineral species, mineral groups, 1. Basic Operation 

plant species, plant groups, man-made objects, etc. The Generally, the first thing that an analyst would do is 

selection is limited only by the extent of spectra placed into display the ternary diagram using raw image data to get an 

storage for use by the computer program. The selection 65 overview of the distribution of the type of ground cover. The 

could be by pointing at a list or typing in one or more names. approximate shape and location of point clusters can be used 

The method of choice is not critical. Both approaches might to infer what the materials are. That might be followed by 
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viplotting the. temaryi diagram from imagcpixels converted to ^- <^iiised for acquiring the. imagery) and plot the point(s) on the . 

either exoatmospheric reflectance or apparent (ground) ternary diagram. Then have the computer program search for 

reflectance to obtain patterns that more closely match the spectrally similar pixels, in the LUT, composed of pixels 

theoretical shapes and locations as shown in FIGS. 1 and 2. converted from either exoatmospheric reflectance or appar- 

The ternary diagram can have class boundaries 5 ent reflectance. It is probable that there would be no exact 

superimposed, optionally displayable, in the manner shown matches. The tolerance can be increased by expanding the 

in FIG. 3 or as implied by FIG. 6, to aid in the visual search to include a region within the ternary diagram 

identificationof the pixel class(es). The interpretation of the enclosed by an ellipse with the major axis parallel to the 

spectral class represented by the plotted points can be ^mineral line' and enclosing the plotted point(s) of interest, 

assisted if the plotted points are colored in relation to their lo The size of the ellipse can be increased repeatedly, 

position, which is governed by their respective ternary Eventually, some target pixels will be identified. They may 

percentages. This can be done by assigning the color be false-positive indications, though. However, it may still 

whereby the red-green-blue components of the resulting hue indicate favorable regions for examination with sensors with 

are directly proportional to the mid-IR, near-IR, and visible higher spatial or spectral resolution, or reconnaissance on 

percentages, respectively. That is, a ternary diagram point 15 the ground. 

with a large TM band 4 percentage will appear greenish, a Alternatively, if one finds it necessary or desirable to work 

point with a large TM band 5 or 7 percentage will appear with the raw image, composed of the DNs as provided by the 

reddish, and a point with a large TM band 1, 2, or 3 data supplier, one can select a point and/or region on the 

percentage will appear bluish. The exact hue will be deter- ground that is of interest. Such a feature might be a known 

mined by the relative proportions of each of the three band 20 mine, known mineral outcrop, known agricultural crop, or 

contributions. any other target of verified character. The selected target will 

The real-time linkage between the displayed image(s) and be displayed on the ternary diagram. Then, by instructing the 

ternary diagram(s) through a look-up table (LUT) allows program to find all points in the image with ternary ratios 

one to specify a point or region in the ternary diagram, and equal to or approximately equal to the original target, new 

identify (with a bright color and/or bUnking pixels) the 25 targets will be identified in the image. As above, an approach 

spatial location(s) in the image that correspond(s) to that would be to place an ellipse, defining a region of interest in 

specified region of interest. This allows one to select a the ternary diagram, around the selected point(s), with the 

particular mineral or type of vegetation and find where in the long axis of the ellipse parallel to the cluster of the target 

image it might occur. In order to account for mineral and class (i.e. vertical for minerals and at 600 for vegetation for 

plant variations in reflectance arising out of compositional 30 this embodiment). The length of the long axis and eccen- 

differences and to compensate for measurement errors, it is iricity can be varied to increase the tolerance of acceptance 

generaUy advisable to plot the desired target material from for display in the image. 

more than one spectral library, if available. The search strategy can be refined further by using more 

The displayed image(s) may be false -color composites of than one ternary diagram at a time. As an example, some 

raw (binary byte DN) data, radiance, exoatmospheric 35 sulfides plot in very different places on a TM band 4,7,3 

reflectance, or corrected apparent-reflectance. However, for combination than they do with a TM band 4,53 combina- 

matching against spectral library entries, the best results will tion. To take advantage of this, one would plot the particular 

be obtained if ternary percentages derived from apparent sulfide of interest from the TM -bandwidth convolved min- 

reflectance or at least exoatmospheric-reflectance pixel val- eral library. A tolerance would be selected by placing a 

ues are loaded into the LUT. 40 bounded area around the point(s) such as previously 

Expanding the target-seeking operation from a single described with the eUipse method above. The computer 

region in the ternary diagram to a plurality of regions, and program will then employ an extended look-up table, incor- 

assigning different colors to the image pixels that corre^ond porating the additional information of the percentages for 

to each of the plurality of ternary diagram regions defined, the second virmal ternary diagram. When an image pixel 

will result in a thematic map. For example, all image pixels 45 falls within the tolerance range of the selected library entry 

that correspond to the vegetation region in the ternary (e.g. sulfide species) for one diagram triad, (e.g. TM 4,73), 

diagram might be colored green, all pixels corresponding to the program will also check to confirm that the image pixel 

minerals might be colored red, pixels intermediate between also is within the tolerance range for the second selected 

the two regions might be colored yellow, and pixels corre- triad (e.g. TM 4,53). Effectively, a logical AND test is being 

sponding to the water region in the ternary diagram might be 50 performed. This approach may be extended to a plurahty of 

colored blue. To achieve finer discrimination, smaller diagram triads. 

regions or segments can be defined, with an inaease in the However, it is not advisable to use most of the possible 

number of display colors used to differentiate the thematic triads for AND-testing for the same reasons that it is not 

classes. This process is analogous to the selection of training advisable to use them singly — principally, the vegetation 

sites in conventional supervised classification, except that 55 and mineral clusters overlap, resulting in poor discrimina- 

the sites arc defined within the ternary diagram feature-space tion and false-positive results. Therefore, the six combina- 

instead of being defined in the spatial-domain feature-space tions shown in FIGS. 1 and 2 are the optimum choices, 

of the image. 3. Mixed-Pixel Problem 

A plurality of hierarchical schema, under computer At Landsat TM resolution, most of the image pixels will 
control, may be devised to automatically classify an image, 60 be mixed. For example, a metamorphic skam is not monom- 
using virtual segmentation boimdaries based on the ternary ineralic; it will be composed of several calc-silicate miner- 
percentages of the optimum band triads; FIG, 6 is one als. An approach to exploration then would be to assume 
example of such a schema. what the two or three most abundant minerals might be, e.g. 
2. Examples of Operation epidote and grossularite and place an ellipse on the Trian- 

One way in which the ternary diagram can be used in 65 gular diagram with the foci of the ellipse coincident with the 

exploration would be to select a particular mineral kom a centroids of epidote and grossularite points, as determined 

spectral library (convolved to the band widths of the sensor from a spectral library for minerals. The image would then 
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be searched . with, variable.:. eccentricity,..ellipseS' for* all ^ .4. Atmospheric Corrections . ^.-^ • •'■^■ii\*.\-^^^^.^~-.-ir^ . » rw - 
matches that fall within the perimeter of the ellipses. This Experience has shown that plotting the ternary percent- 
will allow identification of composition end-members as ages derived from the raw digital numbers (DN) in the 
well as the intermediate composition mixed-pixels. Landsat TM images typically provides useable results. 

Mixed-pixels are always a problem in thematic classifi- 5 However, to achieve improved target identification and more 

cation. Since mineral and vegetation groups are linearized, reliable estimates of proportions, it is desirable to convert 

only groups at the ends of the clusters can be assumed to be the raw DNs to apparent reflectance. This requires using the 

pure. Intermediate points could possibly be the result of calibration metadata supplied with the imagery to convert to 

inhomogeneous reflectors. However, depending on the man- absolute radiance, correcting for sun elevation and distance 

ner in which the groups are ordered, mixing may be min- lO and solar irradiance variation with wavelength, and to finally 

eralogically improbable. Similarly, plant communities com- convert the DNs to normalized exoalmospheric reflectance, 

posed of the vegetation end-members may be improbable for By additionally correcting for atmospheric absorption and 

the particular geographic region, leading to the conclusion scattering, a value approximating actual reflectance may be 

that the pixels are actually representative of the vegetation derived from the exoatmospheric reflectance. Plotting the 

on the ground. is 3 -band ratios derived from the apparent reflectance values 

Ternary Multispectral Analysis may be most accurate provides the most accurate placement of clusters and indi- 

when incorporated with a GIS or an Artificial Intelligence vidual points in the triangular feature -space for comparison 

Expert program such as Prospector, to eliminate the improb- with the idealized locations derived from laboratory-derived 

able interpretations, based on what is known about the area reflectance values. 

and what is known about mineral and plant associations. 20 As with hyperspectral image analysis, accurate correction 

A more serious — and common — ^problem will be mixed- for atmospheric absorption and scattering effects are crucial 

pixels consisting of vegetation and minerals. However, an to accurate and fine discrimination of mineral and plant 

advantage of using Ternary Multispectral Analysis is that species. Ideally, a radiative transfer model (RTM) should be 

mixed-pixels, composed of mineral and plant end-members, used. Lacking ready access to an RTM or the requisite 

are apparent by the intermediate position between the min- 25 environmental data for the time of image acquisition, a 

eral and plant clusters. This effect is most pronounced with first-order correction can be made by correcting for additive- 

the TM band 4, 5, and 3 or 4, 7, and 3 combinations because haze Rayleigh scattering. In the dark-object subtraction 

of the stronglinearization of minerals and distinct separation method described by Chavez (1989), a dark area in the 

of minerals and plants. If either the mineral or the vegetation blue-green or green band — such as a deep, clear lake or a 

end-member can be identified from the presence of pure 30 deeply shadowed gorge — should have a zero (0) radiance 

pixels, or a priori knowledge about the geology or biology, value, or close to it. Subtract the smaUest value in a 

then the other end-member and proportion can be estimated visible -wave length band from all pixel radiance values in 

by extrapolating the line connecting the known end-member that band. Since Rayleigh scattering varies approximately 

and the mixed-pixel. If either end-member is of particular inversely with the 4^ power of the wavelength for clear 

interest, ground reconnaissance of the areas of the mixed 35 skies, the other two bands can be adjusted firom the empiri- 

pixels may be warranted. cally observed bias. Observe the effect on the points plotted 

If one of the mixed-pixel end-members is of no interest, on the ternary diagram. If too large a bias is removed, (fikely 

such as the soil background, and one wishes to derive a with this approach) then the clusters will be shifted away 

thematic map of plant types, the line connecting the mixed- firom the comer with the visible-band ratio on the ternary 

pixel with the local soil end-member can be extrapolated to 40 diagram. In particular, mineral cluster 11 will not fall along 

the plant region on the edge of the ternary diagram and the bisector 12. If no correction is applied (or insufficient bias 

image mixed-pixels can be replaced with pixels correspond- removed) then one can expect the clusters to be shifted 

ing to the DNs of those plants. toward the triangle apex with the visible-band ratio values. 

The ternary diagram also offers a new type of vegetation As shown in FIGS. 1, 2, 3, and 4, that would be the bottom, 

index that is specific to particular plant species or plant 45 right-hand corner of the triangle. An approximate approach 

communities. TTie relative distance of the mixed- pixel along to correcting for Rayleigh scattering would be to correct 

the line connecting the two end members may be used as an only the visible wavelength band used (TM 1, 2, or 3), since 

estimate of the relative proportion of vegetative ground Rayleigh scattering in the IR is usually negfigible. 

cover. The distance between the soil/mineral end-member To empirically establish atmospheric correction 

and the plant end-member is measured. The proportional so parameters, an approach that would be helpful would be to 

distance, along the transect connecting the two end- plot triads that produce linear clusters irrespective of the 

members, of a particular pixel represents the percent of specific surface materials. A plot of TM bands 4, 7, and 5; 

ground cover of the vegetation. If the plant is prone to 2, 7 and 1; and 4, 5, and 3 provide examples of linear clusters 

complete coverage in at least some areas, then pixels will be that can be used for absorption and scattering corrections, 

available to establish the end-member position of the plant. 55 Plots of these triads are shown as FIGS. 7 A, 7B, and 7C. 

If not, then the position can be established by plotting a point Hiesc would serve as guides to indicate if the corrections 

from the spectral library or entered into the program from have been done correctly. If all three bands have been 

field spectrometer measurements. An alternate approach is corrected properly, the point clusters obtained firom the 

to extrapolate any hnear cluster of mixed-pixels and esti- corrected image bands should fall on or within the alignment 

mate the position of the vegetation end-member to be where 60 bands for all three diagrams. 

the extrapolated line intercepts the extreme left edge of An example of how the above can be accomplished would 
vegetation cluster 10. be a software implementation that facilitates the rapid deter- 
The great distance between the pure water point and the mination of atmospheric corrections utilizing a mouse- 
mineral-soil clusters on the ternary diagram provides oppor- cursor-actuated slider on screen. This will adjust gain 
mnity for a sensitive estimate of turbidity and approximate 65 (multiplicative haze) and bias (additive haze) and thus move 
composition of the suspended sediment in a manner similar the plotted clusters in real-time to align the cluster(s) with a 
to above. pre-drawn line(s) or band(s) on the ternary diagram(s) and 
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- then, record the parameteis jequired* to affect the.traDslation 
(s). These parameters will be used in subsequent cluster 
plots for analysis. While bands 1, 2, and 3 are the only ones 
requiring potentially significant corrections for Rayleigh 
scattering, multiplicative haze (Mie scattering) corrections 
for higher bands have proven difficult to deal with. One 
could display three different ternary diagrams (FIGS. 7A, 
7B, and 7C) simultaneously and have them linked together 
in real-time. Thereby, when the slider was adjusted for a 
particular band, one could observe the effect on the plots 
using all three bands of the bands selected for Ternary 
Multispectral Analysis. 
5. Standardized Images 

If, when the ternary percentages are calculated for loading 
43 into the LUT, the values arc also written to a sequential 
file, transformed images will be created. If those percentage- 
images are subsequently displayed in a manner that the 
derivative image obtained from the visible band is directed 
to the blue gun of the CRT (cathode ray tube), the NIR band 
is directed to the green gun, and the Mid-IR band is directed 
to the red gun, one will have created a standardized image 
with a number of desirable properties, such as: 

1) Vegetation will appear in shades of green; 

2) Water will appear in shades of blue; 

3) Most common rock-forming minerals will be gray; 

4) Minerals of economic importance, such as oxides, will 
be reddish; 

5) Borates will be cyan, and; 

6) Copper minerals will be in a magenta region. 
However, the real advantage is that interpreters can 

extract composition information from a false-color image 
that is predictable in color associations. Performing a con- 
trast stretch on the standardized image will increase the 
saturation at the expense of shifting the hue slightly and 
distorting the shape of the ternary diagram clusters. 
Therefore, if the analyst decides to do a contrast stretch, it 
is not recommended that the resulting image be used sub- 
sequently for creating a new ternary diagram. 

CONCLUSIONS, RAMIFICAnONS, AND 
SCOPE 



Accordingly, the reader will see that the Ternary Multi- 
spectral Analysis method and process expands the field of ^5 
multispectral image processing by providing new and dif- 
ferent approaches to identifying targets, creating thematic 
maps, performing atmospheric conrections, unmixing pixels, 
and displaying standardized-color false-color images. 
Highly accurate thematic maps are possible, based on spec- 
tral characteristics of the ground features, using fewer bands 
than typically are employed. 

Furthermore, Ternary Multispectral Analysis has the addi- 
tional advantages that it: 

(a) reduces the level of skill and experience required for 
multispectral analysis; 

(b) allows for automatic hierarchical thematic classifica- 
tion; 

(c) allows interactive generation of thematic maps without 
need for a photointerpreter, ground truth, or more than 
3 TM bands; 

(d) minimizes the impact of the *shadow class' through an 
inherent normalization; 

(e) allows detection of ^ecific spectral targets; 

(f) estabUshes a method to estimate end-member propor- 
tions in a mixed-pixel; 
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■ (g) provides a means to accomplish empirical atmospheric 
corrections; 

(h) provides for displaying false-color images with stan- 
dard colors; 

(i) establi^es a method to calculate species-specific veg- 
etation indices; and 

(j) allows the estimation of relative turbidity, and type of 

sediment in water. 
While my above description contains many specificities, 
these should not be construed as limitations on the scope of 
the invention, but rather as an exemphfication of one pre- 
ferred embodiment thereof, with minor variations. Many 
other embodiments are possible. For example, the use of an 
adaptation of the automatic classification schema (FIG. 6) 
may allow the automatic estimation of cloud cover percent- 
age of an image for quality control purposes. Another 
example is that if an analyst wants to work with a false- 
color, principle components image (Karhunen-Loeve 
transform), notorious for having unpredictable colors, a 
ternary diagram linked to the image will assist the analyst in 
interpreting the type of ground cover. One last example is, 
by spatially linking two ternary diagrams in real-time to an 
image on the computer display, the distinctly mixed-pixels 
in a TM 4, 5, and 3 combination can be high-lighted 
simultaneously on a higher variance, TM 4, 7, and 1 or 2 
combination for estimating proportions of end-members. 

Widely varying embodiments of the invention will sug- 
gest themselves to those having skill in the art without 
departing from the scope and essence of the following 
claims. Accordingly, the scope of the Ternary Multispectral 
Analysis invention should be determined not by the embodi- 
ments illustrated, but by the appended claims and their legal 
equivalents. 
What I claim is: 

1. A method of identifying spectrally-distinctive reflective 
ground targets, represented by a plurality of pixels in mul- 
tispectral imagery, comprising the following steps: 

(a) selecting at least three image bands from a multitude 
of coregistered multispectral bands, wherein said 
selected bands consist of a band in the visible, a band 
in the near-infrared, and a band in the mid-infrared 
portion of the electromagnetic spectrum; 

(b) transformation means for correcting said pixels' origi- 
nal values for sensor calibration, solar spectrum 
characteristics, sensor and illumination geometry, and 
atmospheric scattering and absorption, thereby result- 
ing in values that closely approximate true reflectance, 
herein referred to as apparent reflectance values; 

(c) conversion means for converting plurality of said 
pixels' apparent reflectance values to ternary 
percentages, wherein each member of a triad of per- 
centages is the fractional part, of the sum, of one pixels* 
values from each of said selected bands; 

(d) plotting indicia on a conventional ternary diagram or 
triangular plot, wherein each pixel of the plurality of 
said pixels has an indicium corresponding to said pixel; 

(e) determining a plurahty of spectral-class segment 
boundaries by visual inspection, from a reference ter- 
nary diagram containing plotted points or indicia 
obtained from a spectral library of reflectors of interest 
convolved to the bandwidth of the multispectral sensor 
used, wherein said boundaries are changeable and 
determined by class-resolution requirements of a clas- 
sification task; 

(£) transferring said changeable boimdaries to a triangular- 
shaped, transparent template that is congruent with said 
ternary diagram; and 
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. (g) superimposing said template, over said ternary diagram - - 
to delineate plurality of said spectral-class boundaries; 
whereby, the position of the indicia on the ternary diagram, 
with respect to the plurality of the enclosing boundaries of 
said superimposed template, constrains the identification of 5 
the corresponding image pixels to be members of the 
spectral class plotted on said reference ternary diagram, 
within said enclosing boundaries, thus identifying said pix- 
els. 

2. The method of identifying spectrally-distinctive re flee- 
live ground targets, as set forth in claim 1, step (e) wherein 
the template boundaries determination further includes, 
optionally: 

(a) spectral-class boundary determination means empiri- 
cally derived from actual spectral ground-truth to estab- 
lish the position on a master ternary diagram suitable 
for use with the imagery being analyzed; and 

(b) spectral-class boundary determination means empiri- 
cally derived from photointerpretation-derived ground- 
truth of imagery, to establish the position on a master 
ternary diagram sxii table for use with the imagery being 20 
analyzed. 

3. The general method of identifying spectrally- 
distinctive reflective ground targets, as set forth in claim 1, 
comprising the following software-controlled steps: 

(a) providing an input means to ingest multispectral image 25 
data; 

(b) providing a storage means to store programmatic 
control and data; 

(c) reading said image data into storage; 

(d) providing computer input devices of known type, 30 
thereby enabling an operator to select requisite image 
bands interactively; 

(e) reading a subset of multispectral image data, compris- 
ing said selected bands, into an array; 

(f) reading pixels values from said array; 35 

(g) correcting said pixel values, thereby obtaining appar- 
ent reflectance values; 

(h) displaying said pixel values as a matrix, thereby 
forming an image; 

(i) selecting a region of interest in said image, interac- ^ 
lively; 

(j) converting said apparent reflectance values, within said 

region of interest, to ternary percentages; 
(k) generating a graphical output, on a display device, in 
the form of a ternary diagram containing plotted points 
or indicia corresponding to said pixels; 
(1) displaying pre-calculated, general spectral-class seg- 
ment boxmdaries on said ternary diagram; and 
(m) comparing indicia positions to a table of pre- 50 
determined virtual boundaries to establish finer class 
membership discrimination; 
whereby, the steps recited above result in the semiautomatic 
identification of said selected sub-set of pixels through 
programmatic control of a general-purpose computer. 55 

4. The method of identifying spectrally-distinctive reflec- 
tive ground targets, as set forth in claim 1, wherein step (c) 
said pixel's initial values further include, optionally: 

(a) raw digital numbers, representing a function of radi- 
ance recorded by an overhead sensor. 60 

5. The method of identifying spectrally-distinctive reflec- 
tive ground targets, as set forth in claim 1, wherein step (c) 
said pixel's initial values further include, optionally: 

(a) exoatmospheric reflectance values resulting from a 
processing means to correct for sensor calibration, solar 65 
spectrum characteristics, and sensor and illimiination 
geometry. 



«- 6. The method of identifying spectrally-distinctiverTeflec-w 
tive ground targets, as set forth in claim 1, wherein step (d) 
plotting further includes: 

(a) a coloration means to assign a color to an indicium 
where said color is determined by ternary percentages 
in a manner whereby a red component of a resulting 
hue is directly correlated with one percentage, a green 
component with another, and a blue component with 
the remaining percentage. 

7. The method of identifying spectrally-distinctive reflec- 
tive ground targets, as set forth in claim 1, wherein step (d) 
plotting further includes: 

(a) a coloration means to assign a color to an indicium 
where said color is determined by ternary percentages 
in a manner whereby a red component of a resulting 
hue is directly correlated with a mid-infrared 
percentage, a green component with a near-infrared 
percentage, and a blue component with a visible per- 
centage. 

whereby, said indicium is assigned a hue that approximates 
the natural hue of the ground-cover spectral class the indi- 
cium corresponds to. 

8. The method of identifying spectrally-distinctive reflec- 
tive ground targets, as set forth in claim 1, comprising the 
following software -controlled steps: 

(a) providing an input means to ingest multispectral image 
data; 

(b) providing storage means to store programmatic con- 
trol and data; 

(c) reading said image data into storage; 

(d) providing computer input devices of known type, 
thereby enabling an operator to select the requisite 
image bands interactively; 

(e) reading a subset of multispectral image data, compris- 
ing said selected bands, into an array; 

(f) reading a plurality of pixels values from said array; 

(g) correcting said pixel values, thereby obtaining appar- 
ent reflectance values; 

(h) converting said apparent reflectance values to ternary 
percentages; and 

(i) writing said ternary percentages to a sequential file; 
whereby, a transformed image is created, which is in 
machine-readable form, and can be displayed on a graphics- 
enabled computer; thereby, if the near-infrared percentage 
surrogate is directed to the green gun of a cathode ray tube, 
the mid-infrared percentage surrogate is directed to the red 
gun of a cathode ray tube, and the visible percentage 
surrogate is directed to the blue gun of a cathode ray tube, 
a standardized, false-color image will be displayed that has 
the desirable properties of displaying ground cover in hues 
approximating their natural hues. 

9. A method of locating spectrally-distinctive reflective 
ground targets, represented by a plurality of pixels in mul- 
tispectral imagery, and deriving a thematic map, comprising 
the following software-controlled steps: 

(a) providing an input means to ingest multispectral image 
data; 

(b) providing a storage means to store programmatic 
control and data; 

(c) reading said image data into storage; 

(d) providing computer input devices of known type, 
thereby enabling an operator to select requisite image 
bands interactively; 

(e) reading a subset of multispectral image data, compris- 
ing said selected bands, into an array; 



06/16/2004, EAST Version: 1.4.1 



us 6,356,646 Bl 



23 



24 



30 



; (f) reading, pixels ^values fromc said^array;. ^-^ — . ^ . 

(g) correcting said pixel values, thereby obtaining appar- 
ent reflectance values; 

(h) displaying said pixel values as a matrix, thereby 
forming an image; ^ 

(i) converting said apparent reflectance values to ternary 
percentages; 

(j) loading a lookup table, sequentially, with said ternary 
percentages; 

(k) generating a graphical output, on a display device, in 
the form of a ternary diagram containing plotted points 
or indicia; 

(1) segmenting said ternary diagram by selecting a plu- 
rality of boundaries, interactively, that defines a desired is 
spectral class; 

(m) comparing entries in said lookup table with bound- 
aries of spectral-class segment; 

(n) highlighting image pixels if said pixels correspond to 
said lookup table entries within said region of interest; 20 
whereby, the steps recited above result in locating the 
plurality of said pixels in said image that spectrally corre- 
spond to said indicia in said selected spectral-class segment; 
further, the steps (1)> (tn) and (n), recited above, iteratively 
performed on the whole of the pluraHty of indicia, with 25 
non-overlapping segments, result in a plurality of conespon- 
dences between said segments and the plurality of said 
image pixels; 

whereupon, assigning unique a color to each spectral class, 
for which there is at least one correspondence in the image 
pixels, results in the generation of a thematic map, depicting 
the types of ground cover as represented by the plurality of 
said spectral classes. 

10. A method of locating spectrally-distinctive reflective 
ground targets and deriving a thematic map as set forth in 
claim 9, wherein steps (h), (J), (k), (m), (n), and (o) are 
omitted and (1) is replaced with the further steps: 

(a) creating a pre-determined plurality of virtual spectral- 
class boundaries, comprising a plurality of test 
parameters, each parameter consisting of no more than 
two ternary percentages, which establish a boundary 
between two super-classes; 

(b) performing hierarchical logical tests between the 
relative values of said virmal, spectral-class boundaries 
and the corresponding ternary percentages for the plu- 
rality of said image pixels; 

(c) continuing said logical comparisons, with additional 
parameters, on the pixel along a decision tree until an 
identification is made; 

(d) assigning an identification label, such as a number, to 
the pixel; 

(e) writing said label to a sequential file; 

(f) iterating the procedure by selecting the next pixel in 
sequence; 

(g) continuing until all of the plurahty of said pixels have 
been tested; 

whereby, the steps recited above result in an automatic 
classification of said multispectral image suitable for 
unattended, batch processing in a volume production 
envirotmient, wherein the resulting thematic map is in a 
machine-readable form. 

11. A method of locating spectrally-distinctive reflective 
ground targets as set forth in claim 9, wherein step (1) is 
replaced with the further steps: 

(a) storing a reference spectral library convolved to appro- 
priate baodpasses; 
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-- iK(b) providing an interactive selection means for extracting.^- 
at least one member from said reference spectral 
library; 

(c) converting said selected spectral hbrary member 
reflectance values to ternary percentages; 

(d) providing a means for establishing a tolerance of 
acceptance of a match; 

(e) comparing said spectral library member, 
consecutively, with every entry in said lookup table; 

(f) highlighting the corresponding image pixel if said 
match is obtained; 

(g) repeating step (e) for each additional said member of 
said spectral library selected; 

whereby, the steps recited above result in locating the 
plurality of said pixels in said image that spectrally corre- 
spond to said selected member of said reference spectral 
library. 

12. A method of interactively, empirically correcting 
multispectral imagery for multiplicative and additive haze 
and atmospheric absorption, comprising the following 
software-controlled steps: 

(a) providing an input means to ingest multispectral image 
data; 

(b) providing storage means to store programmatic con- 
trol and data; 

(c) reading said multispectral image data into storage; 

(d) providing computer input devices of known type, 
thereby enabling an operator to select the requisite 
image bands interactively; 

(e) reading a subset of multispectral image data, compris- 
ing said selected bands, into an array; 

(f) reading pixels values from said array; 

(g) converting raw image pixel-values to ternary percent- 
ages; 

(h) loading a lookup table, sequentially, with said ternary 
percentages; 

(i) duplicating said ternary percentages in a shadow 
lookup table; 

(j) generating a graphical output, on a display device, of 

a plurality of different ternary diagrams containing 

plotted points or indicia; 
(k) displaying at least one pre-determined alignment band 

on each of said plurality of ternary diagrams; 
(1) providing interactive input means whereby an operator 

adjusts visible indicators of gain and bias; 
(m) providing a means to record and apply adjusted gain 

and bias values to said ternary percentages; 
(n) providing a means to update said shadow lookup table 

with the adjusted percentages; 
(o) regenerating said graphical output comprising 

adjusted said indicia, on said display device, using the 

updated shadow lookup table; 
(p) continuing said interactive adjustments until satisfac- 
tory alignment with plurality of said alignment bands is 

achieved; 

(q) providing a means to convert final, adjusted ternary- 
percentage gain and bias values to values to apply 
directly to the raw data to obtain apparent reflectance; 
whereby, said raw image data is corrected for a multitude of 
multiplicative and additive effects that arc essential to take 
into account to create apparent reflectance images. 
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